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We realize a gapless Majorana Orbital Liquid (MOL) using orbital degrees of freedom and also an SU(2)- 
invariant Majorana Spin Liquid (MSL) using both spin and orbital degrees of freedom in Kitaev-type models 
on a 2-leg ladder. The models are exactly solvable by Kitaev's parton approach, and we obtain long-wavelength 
descriptions for both Majorana liquids. The MOL has one gapless mode and power law correlations in energy 
at incommensuare wavevectors, while the SU(2) MSL has three gapless modes and power law correlations in 
spin, spin-nematic, and local energy observables. We study the stability of such states to perturbations away 
from the exactly solvable points. We find that both MOL and MSL can be stable against allowed short-range 
parton interactions. We also argue that both states persist upon allowing Z2 gauge field fluctuations, in that the 
number of gapless modes is retained, although with an expanded set of contributions to observables compared 
to the free parton mean field. 



I. INTRODUCTION 

Since the experimental realization of gapless quantum spin 
liquids (QSL)ii— in two-dimensional (2D) organic com- 
pounds K-(ET)2Cu2(CN)3 and EtMe3Sb[Pd(dmit)2]2ri— 
there have been many theoretical proposals— -24 for such in- 
triguing phases. Among them, the proposal of an SU(2)- 
invariant Majorana spin liquid (MSL) by Biswas et ai— is 
fascinating and in need of more careful consideration. 

In an earlier work, we constructed an exactly solvable mi- 
croscopic model^ in Kitaev's spirit''^ to study the properties 
of such SU(2)-invariant MSL with Fermi surfaces of partons. 
However, we allowed very low symmetries — lack of parity, 
inversion, and time reversal symmetry (TRS) — to sidestep 
discussing possible perturbations such as Cooper pairing in- 
stability, which can destabilize the gapless QSL phases away 
from the exactly solvable limit. In order to study the sta- 
bility of such new class of gapless QSL and further explore 
their properties, we realize such states on a 2-leg square lad- 
der and show that they represent new quasi-one-dimensional 
(ID) phases. 

We first consider a gapless Majorana Orbital Liquid (MOL) 
realized in a Kitaev-type model on the 2-leg ladder using 
orbital degrees of freedom. The system can be reduced to 
one species of Majorana fermions coupled to background Z2 
gauge fields such that it is exactly solvable and has gapless 
partons with incommensurate Fermi wave vectors. We formu- 
late a long-wavelength description in terms of right-moving 
and left-moving complex fermions Jh/l and show that local 
energy observable has power law correlations at incommen- 
surate "2fci?" wavevectors. Going away from the exactly solv- 
able point, we first consider allowed residual parton interac- 
tions and find that there is only one valid four-fermion term 
and it is strictly marginal; hence, the MOL is stable to such 
perturbations. 

An important question is the stability of the MOL to allow- 
ing Z2 gauge field fluctuations, as these lead to confinement 
of partons in gapped phases in so-called even Z2 gauge the- 
ories in (lH-l)Dfiii2Z We argue that because of the nontrivial 
momenta caried by the gapless partons, there is a destruc- 
tive interference for Z2 vortices (instantons) in space-time, 
and hence these are suppressed and do not affect the count 



of gapless modes. The local energy observables obtain new 
contributions beyond the mean field, and in this sense the par- 
tons become "less free", but their bosonized fields still remain 
very convenient for characterizing the MOL phase. 



We next realize an SU(2)-invariant Majorana Spin Liquid 
(MSL) using both spin- 1/2 and orbital degrees of freedom^^ 
at each site of the 2-leg ladder4£ The system can be reduced 
to three species of Majorana fermions coupled to background 
Z2 gauge fields such that it is exactly solvable and has gapless 
partons with incommensurate wave vectors. We formulate 
long-wavelength description in terms of three right-moving 
and left-moving complex fermions (.f^/^.' /fl/L' ^r/l) 
transform as a vector under spin rotation. Because there 
is no global U(l) symmetry, in addition to familiar four- 



fermion residual interactions expressed as f^^ fj^^ f]{ fl, there 
are other allowed terms such as ff^ fi} f]i f^j^ ■ Despite of 
having more allowed interactions, a weak coupling renormal- 
ization group (RG) analysis gives a large regime of a stable 
phase. Similarly to the MOL case, we argue that such MSL 
with gapless matter can be also stable against Z2 gauge field 
fluctuations even in (iH-l)Di^iliii^ 



The paper is organized as follows. In Sec. [Ill we realize 
the MOL with one fermion species in a Kitaev-type model^ 
on the 2-leg ladder and consider its long-wavelength proper- 
ties and stability against perturbations. In Sec.|IIIl we real- 
ize the SU(2) MSL and use weak coupling RG analysis to 
study the stability of such phase against residual parton in- 
teractions and also discuss the stability against gauge field 
fluctuations. We conclude in Sec. |IV] with some discussions. 
In Appendix |A] we consider more abstractly the stability of 
gapless U(l) matter against Z2 gauge field fluctuations in 
(1h-1)D. In AppendixiBl we give long-wavelength description 
of the SU(2) MSL and discuss observable properties. In Ap- 
pendix |Cl we consider Zeeman magnetic fields on the SU(2) 
MSL. In Appendix ini we realize the SU(2) MSL in a model 
with explicitly broken time reversal symmetry and show that 
this case has a larger window of stability to weak perturba- 
tions. 
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II. GAPLESS MAJORANA ORBITAL LIQUID (MOL) ON A 
TWO-LEG LADDER 



We begin with a "spinless" (one species) MOL realized in 
a Kitaev-type model on a 2-leg ladder shown in Fig. |l(a)| The 
Hamiltonian is 



>unit cell 



where 



Ho — 



A-link,(ifc) 



'''l ■'"2 "^3 "^4 I 



l^n^, = T2 Ti T4 Tg . 



(1) 



(2) 

(3) 
(4) 



The r Pauli matrices can be thought of as acting on two-level 
orbital states. The Wp terms, with p = U^z or Dj^j, formed by 
X and z ox y and z links respectively, are plaquette operators 
which commute among themselves and with all other terms in 
the Hamiltonian and are added to stabilize particular flux sec- 
tor, see Fig. |l(a)| Following Kitaev's approach, we introduce 
Majorana representation as 



(5) 



with the constraint = b^b^bjCj — 1. The Hamiltonian can 
be rephrased as 



Ho = i'Y'^okJjkCjCk , 

{]k)ep 



(6) 
(7) 



where u,- 



ijk = —ibjb'^ for A-link (jk) and the product in the 
last line is circling the plaquette. 

Following familiar analysis in Kitaev-type models, we ob- 
serve that in the enlarged Hilbert space, ujk commute among 
themselves and with the Hamiltonian, and we can proceed by 
replacing them by their eigenvalues ±1 and interpreting as 
static Z2 gauge fields. The Wp terms, with Kp > assumed 
to be sufficiently large, can be used to stabilize the sector with 
zero fluxes through all elementary plackets, and this can give 
a gapless phase. In our work, we fix the gauge by taking 
Ujk — 1 for bonds j ^ k as shown by the arrows in Fig. |l(a)| 
There are four physical sites per unit cell, so there are four 
Majoranas per unit cell. From now on, we replace the site 
labeling j with j = {X, a}, where X runs over the one- 
dimensional lattice of unit cells of the ladder and a runs over 
the four sites in the unit cell, see Fig. |l(a)| The Hamiltonian 
can be written as, 

H ^'Y^j-^jkCk = Y, Cx,aAx,a;X',a'CX',a'- 

i'jk) {(X,a),{X',a')) 

There is translational symmetry between different unit cells, 

and Ax,a,X',a' = Aaa'{X - X'). 

In order to give a concise long-wavelength description, it 
will be convenient to use familiar complex fermion fields. To 
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FIG. 1. (a) Graphical representation of the exactly solvable Kitaev- 
type model on the 2-leg ladder and its solution in the zero flux sec- 
tor. The c Majoranas propagate with pure imaginary hopping am- 
plitudes specified by the couplings Jx, Jy, Jz, and J',; the signs 
in our chosen gauge are indicated by the arrows and the four-site 
unit cell is also indicated, (b) Dispersion of complex fermions that 
solve the Majorana problem for parameters {J^, Jy, Jz, J'z} = 
{1.2, 0.8, 1.0, 1.1}. 



this end, we can proceed as follows. For a general Majorana 
problem specified by an anti-symmetric pure imaginary ma- 
trix Ajk, we diagonalize Ajk for spectra, but only half of 
the bands are needed while the rest of the bands can be ob- 
tained by a specific relation and are redundant. Explicitly, for 
a system with 2m bands, we can divide them into two groups. 
The first group contains bands from 1 to m with eigenvector- 
eigenenergy pairs {v},^k,^b,k}, where b = l,2,...,m are 
band indices, and the second group contains bands from m + 1 
to 2m related to the first group, {vbi=rn+b,k,^b'=ra+b.k} = 
{^1 _fcj — eb,-fe}- Using only the bands with b = 1 to to, we 
can write the original Majoranas in terms of usual complex 
fermions as 



c(X, a) 



2 



E E [e''''«6,fe(a)/fc(fc)+H.c.], 



fc=l keB.Z. 



where Nuc is the number of unit cells, B.Z. stands for the 
Brillouin Zone, and the complex fermion field / satisfies the 
usual anticommutation relation, {fh{k),fb'{k')} = 5bb'6kk'- 
In terms of the complex fermion fields, the Hamiltonian be- 
comes 



H = E E 2£''(^) flik)fb{k)-- 



6=1 fceB.Z. 



(8) 



In the present case, 2m = 4 and therefore two bands are suf- 
ficient to give us the full solution of the Majorana problem. 

The above approach can be applied to any general Majorana 
problem and is needed when we consider a model lacking any 
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symmetries in AppendixiDl In the present case, we require the 
model to respect time reversal symmetry^ and leg interchange 
symmetry, which allows us to introduce convenient complex 
fermion fields already on the lattice scale as follows 



MX) = 

fuiX) = 



c{X,l) + ic{X,4) 
2 

-ic{X,2) + c{X,3) 



(9) 
(10) 



The Hamiltonian becomes 



X ^ 

- [j,fl{X)fu{X) + Jyfl(X)h{X + 1) + 



H.c. 



where we ignored constant contribution. It is easy to calculate 
the band dispersions. 



e(fc) = J+ ± ^ ( J7)2 + J2 + J2 + 2J^Jy cos(fc) 



(11) 



with J: 



± _ 



[Jz ± J'z)!"^- The spectrum is gapless for 
\Jx — Jy\ \J JzJ'z < Jx + Jy, whcrc without loss of gen- 
erality we assumed all couplings to be positive. For an illus- 
tration of the energy spectrum, we take { , Jy, J^, J'^} = 
{1.2, 0.8, 1.0, 1.1} and show the two bands of the complex 
fermions in Fig. |l(b)| labeled from top to bottom as band-1 
and band-2. We note that the gapless phase occur in a large 
parameter regime and there is no fine tuning here. The specific 
parameters are chosen to emphasize that we do not require any 
symmetries other than time reversal and leg interchange. 

The band-2 crosses zero at kpR and kpL = —kpR from 
time reversal. For long wavelength physics, we can focus 
on this band and introduce continuum complex fermion fields 
fB./L \ for the lattice Majoranas, we obtain the expansion. 



c{X,a)^ b 

P^R/L 



«2.p(a)/p(X)+H.c.] .(12) 



From the detailed band calculation, at the right Fermi point 



4J^ 



(13) 



where ^ = (J^. + Jyc'^^'''^) / J'^. Using time reversal 
invariance, for the left Fermi point we get V2,l{o) ~ 
^(a). The effective low energy Hamiltonian den- 
sity is 



n = VF 



(14) 



describing a one-dimensional Dirac particle with Fermi ve- 
locity vf = JxJy sm{kpB^) / Jt ■ We list the symmetry trans- 
formations of the continuum fields in Table |T] (ignoring the 
"spin" indices there). In particular, the leg interchange sym- 
metry prohibits terms of the form /r/l from the continnum 
Hamiltonian that would gap out the spectrum. 



Fixed-point tlieory of Majorana orbital liquid and 
observables 



In this subsection, we first give the fixed-point theory of 
the MOL and then we will consider bond energy operators 
to characterize such gapless phase. We use Bosonization, 
re-expressing the low-energy fermion operators with Bosonic 
fields42.^ 



with canonical conjugate boson fields: 

[(p{x),0{x')] =iTTe{x~'X') , 



(15) 



(16) 
(17) 



where 9(x) is the Heaviside step function. 

The fixed-point bosonized Lagrangian of such gapless 
MOL is 



-MOL 



1 

2^g 



(18) 



For free fermions, g = 1 and v ~ vp, the bare Fermi 
velocity. Later when we discuss the stability of such a 
phase in Sec. IIIBI we will see that there is only one strictly 
marginal interaction which introduces one Luttinger parame- 
ter g. To detect the gaplessness of the phase using physical 
(gauge-invariant) observables, here we consider bond-energy 
operatorsjiS i5*/°(X), which we further categorize into sym- 
metric or anti-symmetric with respect to the leg interchange 
symmetry. The specific microscopic operators are 

B'/'^iX) = r^(X, l)r^(X, 2) ± r^(X, 4)t^(X, 3) 

= iui2c{X, l)c{X, 2) ± iu43c{X, 4)c(X, 3), (19) 

where we used Majorana representation, Eq. (|5]l. In our gauge, 
after expansion in terms of the continuum complex fermions 
using Eq. ( fT2b . the Fourier components are organized as fol- 
lows 



Q=0 



iIjr + iUl = — , 



Bl 



Bt 



kFR-\-kFL 



fl.fR = 
fhfR = -ie 



i2lp 



(20) 

(21) 
(22) 



[Note that with TRS, the wave vector kpp + /cfl is the same 
as Q = 0; to be more precise, we should write a Hermitian and 
time reversal symmetric combination, Sq^q = «/l/_r + H.c.] 
Thus, the symmetric bond-energy correlations are expected 
to decay with oscillations at incommensurate wave vectors 
±2kpii, while the anti-symmetric bond-energy correlations 
decay without oscillations. Such a sharp difference can be 
confirmed in exact numerical calculations. 

In the bosonized form, the scaling dimension of each term 
is apparent. 



A[S^.o] = 1, 



(23) 
(24) 

(25) 
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|<fi'(0)S'(X)>| 




(a)Leg-symmetric energy correlation 

\<ByO)B''iX)>\ 




(b)Leg-anti-symmetric energy congelation 

FIG. 2. Figures (a) and (b) illustrate power law behaviors of the sym- 
metric and anti-symmetric bond energy correlations, with B*/" de- 
fined in Eq. l |19t . in the exactly solvable model with non-interacting 
partons. The system has 500 unit cells and we use the same param- 
eters as in Fig. [T] We plot absolute values and indicate the sign with 
filled circles (blue) for positive correlations and open square boxes 
(red) for negative correlations. The log-log plots clearly show X^^ 
decay (straight lines) with incommensurate oscillations in the sym- 
metric case and no oscillations in the anti-symmetric case. The char- 
acteristic wavevectors can be determined from the structure factor 
study shown in Fig. [3] 



In the non-interacting parton limit, g ^ 1, we expect to see 
all components of bond-energy correlations decay as X^^. 

For illustration, we calculate correlations in the exactly 
solvable model, taking the same parameters as in Fig.[T] Fig- 
ure [2(a)] shows log-log plot of symmetric bond-energy corre- 
lations in a finite system with 500 unit cells, while Figure [2(b)| 
shows anti-symmetric bond-energy correlations."*^ We can see 
the overall X^^ envelope in both figures and also incommen- 
surate oscillations in the symmetric bond-energy correlations, 
which confirm the theoretical analysis above. 

Power-law correlations in real space correspond to singu- 
larities in momentum space, which we can study by consider- 
ing the corresponding structure factors. Figure [3(a) shows the 
symmetric bond-energy structure factor and Fig. 3(b)] shows 
the anti-symmetric bond-energy structure factor It is clear 
that the singularities in the symmetric case occur exactly at 
Q = and Q = ±{kFR — kpL) = ±2/ci?ii = ±2kF (which 
we also mark using values obtained by extracting the Fermi 
points of band 2), while there is only Q = singularity for 
the anti-symmetric case. 

Let us now consider some other operators similar to generic 
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(b)Leg-anti-symmetric energy structure factor 

FIG. 3. Figures (a) and (b) illustrate the symmetric bond-energy 
and anti-symmetric bond-energy structure factors corresponding to 
Figs. |2(a)[ and |2(b)] respectively. Both cases clearly show a singular- 
ity at Q = 0, while the symmetric case also shows singularities at 

±2kF. 



XYZ energy terms but not present in the exactly solvable 
model; this will be also useful for the subsequent discussion of 
the MOL stability. First, operators like r2'(X, l)Ty{X, 2) and 
T^{X, 1)t^{X, 2) have ultra-short-ranged correlations as they 
contain unpaired localized &-fermions. It is more interesting 
to consider operators like {X, 1)t^ {X, 4) defined on the z- 
type (vertical) links in Fig. [T| In this case, even though the 
local operator contains unpaired 6-Majoranas, in the physical 
Hilbert space these can actually be paired at the expense of 
introducing a string product of the gapless c-Majoranas. For 
example, consider calculating correlation between rungs at X 
and X': 

^■(X,^') = t^(X,1)t^(X,4) t^(X',1)t^(X',4) = 
= H (-l)c(X",l)c(X",4)c(X",2)c(X",3)x 



X<X"<X' 

n 

((X,1),(X,4)) < A-link (tj) < {{X'S),{X'A)) 



ut, (26) 



where the last product contains all links on the ladder located 
between the two vertical links excluding {{X, 1), {X, 4)) and 
including {{X' , 1), {X', 4)) and oriented as shown in Fig. [T] 
The second line is 1 in our chosen gauge, and we then have a 



|F(X)| 

0.3' 



F(q) 




FIG. 4. Figure illustrates power-law behavior of the correlation 
T(X~X') = (jr(X,X')),definedinEq.(|26ll. The system has 100 
unit cells in chain length and the same parameters as in Fig.|T] We 
show the absolute values of and indicate the sign with filled 

circles (blue) for positive correlations and open square boxes (red) 
for negative correlations. The log-log plot clearly shows X~^^^ en- 
velope (straight line in the figure). The irregular behavior is due to 
incommensurate oscillations. 



factor of 



i-l)c{X",l)ciX"A)c{X",2)c{X",3) = 

= e^^lfHx")fiiX") + fiiiX")fni^")] 



(27) 



for each unit cell, where we used Eqs. (l9]l-(fT0l). In the present 
gauge, we can write schematically {X , 1)t^ {X , 4) ^ 
nx"<x(-l)c(^"-l)c(^",4)c(X",2)c(X",3), and see 
that this contains non-Iocal Jordan-Wigner-Iike string operator 
in terms of the gapless partons. In the bosonization language, 
the string operator becomes 



n ' 

X"<X 



„i7r[/t(X")/i(X")+/i\(X")/ii(X")] ^ ±z[e{X)+7rnX] 



(28) 



This has scaling dimension 1/4 in the free-fermion case and 
hence the above correlation decays as X~^^^ power law and 
oscillates at wavevector irfi = kp from Fig. |l(b)| It may 
seem unusual that this appears to contain the specific gauge- 
dependent quantity kp', note, however, that in the full calcula- 
tion we used the specific gauge to set the last line in Eq. ( |26] | 
to unity, and the final result is independent of the gauge. 

Evaluating expectation value of the string operator in the 
free fermion ground state leads to a Pfaffian of a matrix 
formed by the Majorana contractions and can be easily com- 
puted numerically for reasonable sizes."*** The results are 
shown in Fig. |4] for a system with 100 unit cells4i The cor- 
responding structure factor is shown in Fig.|5] We can clearly 
see the singularities at ±kp and confirm our theoretical anal- 
ysis. 



B. Stability of Majorana orbital liquid 

Let us now consider going away from the exactly solvable 
point. First, we consider perturbations that are local in the 
continuum fermion fields. This ignores fluctuations in the Z2 
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FIG. 5. Structure factor corresponding to Fig.|4l we also mark the 
expected locations of the singularities, ztkp- 



gauge fields, and we will address stability against confinement 
shortly. In the language of usual complex fermions, there is 
only one vaUd 4-fermion interaction. 



H 



int 



ufifRflfL 



(29) 




This interaction is strictly marginal, and therefore the gapless 
MOL is stable also with 'Hint and has one gapless mode. This 
interaction will renormalize the Luttinger parameter and the 
Fermi velocity to be 



(30) 
(31) 



which completes our description of the fixed-point theory in 
Eq. ( fTSl l and will modify the power laws of various correla- 
tions as discussed above in Sec. Ill Al 

We now want to address the issue of confinement, more 
precisely, the stability of the MOL theory when we allow 
fluctuations in the Z2 gauge fields. As we discuss in Ap- 
pendixlAl allowing Z2 gauge field fluctuations in the (1h-1)D 
space-time is like allowing half-vortices in the phase field in 
the bosonized harmonic liquid description and corresponds to 
allowing terms Ai/2Cos(6' + kpX + ai/2) in the dual har- 
monic liquid description, Eq. ( fTSl l. The key point is that this 
term is oscillating for generic kp and hence averages out to 
zero (the underlying physics is destructive interference due to 
Berry phases). Thus, our gapless MOL with incommensurate 
momenta carried by the fermion fields persists also in the pres- 
ence of Z2 gauge field dynamics even in (1h-1)D, in the sense 
that we retain the gapless mode. 

One may worry about the precise connection between the 
present system and the schematic Z2 gauge theory plus U{\) 
matter at incommensurate density considered in AppendixlAl 
Indeed, the connection is only crude, and we do not have one- 
to-one correspondences. Nevertheless, we can bolster our ar- 
gument by considering explicitly some allowed perturbations 
to the exactly solvable model. Consider, e.g., adding small 



6 



general XYZ interactions Y.{^J)Y.^.=x,y,.J'Jij'ri'^j on all 
bonds in a manner respecting the underlying lattice symme- 
tries. As we have discussed earlier, SJ^'^^ terms on the x-type 
bonds and 5 J^'^ terms on the 2;-type bonds have short-range 
correlations and hence constitute irrelevant perturbations (of 
course, they can renormalize the Luttinger parameter). On the 
other hand, SJ^ ^ terms on the z-type bonds have power law 
correlations. However, these correlations oscillate at the in- 
commensurate wavevector, see Fig. |4] and Fig.|5] Hence such 
terms, whose structure is similar to A1/2 cos(6'+fci?X+ai/2), 
cf. Eq. ( |28] |, are washed out from the low energy Hamiltonian. 
Thus, the fixed point description is the same as described ear- 
lier, but with the additional remark that now generic energy 
correlations that are symmetric under the leg interchange will 
also obtain a contribution oscillating at wavevector kp with 
scaling dimension g/A. 

Finally, we remark that the Z2 gauge fluctuations do lead to 
confinement in our 2-leg model in gapped regimes, e.g., when 
the Jz terms dominate over the J^, Jy terms in the original 
Hamiltonian Eq. ([T]i. In this regime, we can start with effec- 
tive (super)-spins on the rungs formed by the large terms 
(e.g., after conveniently making the J, coupling ferromag- 
netic). We perturbatively derive effective Hamiltonian gov- 
erning these effective spins, which works out to be an Ising- 
like chain and has two degenerate ground states. Adding the 
Sj^^y perturbations on the z-type bonds gives local longitudi- 
nal fields in this Ising chain and immediately lifts the degen- 
eracy. Hence, there is a unique ground state. 

Furthermore, creating a single domain-wall-like excitation, 
which behaves as a free particle in the exactly solvable model, 
requires infinite energy in the presence of the longitudinal 
field. On the other hand, a pair of domain walls, kink and anti- 
kink, are allowed, but to separate one from the other requires 
energy linearly proportional to the distance between them. 
Therefore, such SJ^'^ perturbations on the 2-type bonds give 
linear confinement of particles that were free at the exactly 
solvable point, and this applies to all particles that carry gauge 
charge with respect to the Z2 gauge field in the exactly solv- 
able model. 



in. GAPLESS SU(2)-INVARIANT MAJORANA SPIN 
LIQUID (MSL) ON THE TWO-LEG LADDER 

We now want to consider Majorana spin liquids with more 
degrees of freedom, in particular with physical spin degrees 
of freedom, and see what new issues and features arise in this 
case. In order to construct spi n SU(2)-invariant Kitaev-type 
model, we follow Refs. I35ll38l and l39l to take a system with 
both spin and orbital degrees of freedom on each site. The 
complete Hamiltonian is 

where 

A-link,(jfc> 



H'q is a Kugel-Khomskii-like Hamiltonian with it being the 
spin- 1/2 Pauli matrices and f being the Pauli matrices act- 
ing on the orbital states, while the Wa^. and Wa^^ terms are 
given in Eqs. (O-dHi. 

Introducing Majorana representation of spin- 1/2, we write 
the spin and orbital operators as 

T- = -^^e"^-<dJ. (35) 

On each site j of the 2-leg ladder, we realize the physical 
four-dimensional Hilbert space using six Majorana fermions 
c^, c'a, rf, (Fa, (K, and df, with the constraint D,- = 
—iCjC^Cjd^djdj = 1 (namely, for any physical state 
|*)phys, we require Dj|$)phys = |3')phys)- Therefore, 
a"T^\^)phys = «c"d^|$)phys- In terms of the Majoranas, 
the Hamiltonian can be rephrased as 

n'o = tJ2ujkJjk (36) 

(jk) a=x,y,z 

and the Wp terms are the same as in Eq. (|7|l with iijk = 
-id^d^ for A-link (jk). 

For long-wavelength description, much of the development 
in Sec.|II]can be directly applied here with the replacement, 
c — > c", / — > a = x,y, z. We now have three fermion 
species with identical dispersion taken to be similar to that 
in Fig. |l(b)[ and we introduce right and left moving complex 
fermion fields j'^jj^ as in the spinless case. Under SU(2) spin 
rotations, the triple ^^■y-'^ transforms in the same way as the 
physical spin a^-"^'^ . 

Just as in the MOL case in Sec|IIl we first estabUsh the fixed 
point structure ignoring the gauge field fluctuations. In or- 
der to study the stability of such gapless SU(2)-invariant Ma- 
jorana spin liquid under weak perturbations, we write down 
most general four-fermion interactions and perform Renor- 
malization Group (RG) studies. The allowed four-fermion in- 
teractions are highly constrained by symmetry. In addition 
to the SU(2) spin rotation invariance, these terms must be 
preserved by Projective Symmetry Group (PSG)^ of spatial 
translational symmetry, time reversal symmetry, and leg in- 
terchange symmetry. We list the symmetry transformations 
in Table H] and write the allowed non-chiral interactions (i.e., 
connecting right and left movers) as 

Hint = UpJbJl - U„iJr ■ Jl + Ua2lliL^RL 

+ WiilRLlRL + il-C.), (37) 

where we defined 

Jp = J2fp-fp^ (38) 

a 

j7i? = -z^e"^Vf/^ (39) 
lRL=J2fRfL- (40) 
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TABLE I. PSG transformation properties of the continuum fields 
under Tx (spatial translation symmetry), Q (time reversal transfor- 
mation plus gauge transformation);^ AI (leg interchange transfor- 
mation plus gauge transformation). We also note that under spin 
rotation, /p = {ff., fp, fl.) and fl = {fp fp, fp) trans- 
form as 3-dimensional vectors. Note that below, P = R/ L and 
P = -P ^ L/R 







e 


M 


fp ^ 




fp; i 


-^fp 






fp , -i 





The general expression 'Hint in Eq. ( |37] ) contains familiar- 
looking four-fermion terms II' fiifL '^hat conserve 
fermion number, and also terms f]^.fpf2.fL ^^^^ 
serve the fermion number but are nevertheless allowed by all 
symmetries of the problem. The less familiar terms need to be 
considered since the microscopic Majorana Hamiltonian does 
not have U{1) particle conservation, which is a new feature in 
such Majorana liquids. 

We remark that the time reversal and translation symmetries 
alone would allow yet other terms expressed as /aVa/^/l 
and in fact would also allow a bilinear term {Hbl + H.c.) 
in the Hamiltonian that would immediately open a gap in the 
spectrum. However, these terms are prohibited if we also re- 
quire the leg interchange symmetry, which is hence crucial for 
the time-reversal invariant SU(2) MSL. 

The weak-coupling differential RG equations are 

Up = (m^2 + 2ucriWo-2 - 4?i;|) , (41) 

ilal = 7; (-"ctI + 2UaiUa2) , (42) 

Ittv 

U(j2 = -Z (-3u^2 ^ GUalUcr2 - 4:wf) , (43) 

ZlTV 

W4 = TT— (^2it^i - 4it^2 - 4up) u;4, (44) 

where v is the Fermi velocity of right and left movers and 
O = dO/dl with £ being logarithm of the length scale. The 
only fixed points have u*„i = w*2 = W4 = 0. Stability to 
small deviations in u'4 requires w* > 0. If we consider small 
deviations in u^i and Ucr2 setting wX — 0, the RG equations 
can be written as 

9p = 3Mp + u^2 = 0, (45) 

Ual = + 2lio-lWcr2), (46) 

ZTTW 

Ua2^-Z (-3w^2 ~ 6Wo-lU^2), (47) 

and the last two equations are essentially identical to the RG 
equations in a level-one SU(3) Wess-Zumino-Witten (WZW) 
model discussed by Itoi and Kato.=52 Translated from their 
analysis, the stability to small deviation in Ud and u„2 re- 
quires u„i > 0, Uo-i + Ucr2 > 0. In a stable flow, Up reaches 
some fixed value, w* > 0, and is strictly marginal; Ug-i and 
Ucr2 approach zero from the specific region described above 



and are marginally irrelevant; finally, u'4 flows to zero as long 
as M* > and is irrelevant. Thus, we have one Luttinger pa- 
rameter in the "charge" sector In Appendix iBl we give the 
fixed-point theory of the SU(2) MSL and list observables that 
can be obtained as fermion bilinears. We find that spin op- 
erator, Eq. ( IB8I 1, spin-nematic operator, Eq. dBlOt , and bond- 
energy operator, Eq. ( |B9t . have correlations that decay in a 
power law with oscillations at incommensurate wave vectors, 
which is one of the hallmarks of such Majorana spin liquids 
as we discussed in Ref . [35l in a 2d example. 

The inclusion of the Z2 gauge field fluctuations in this 
quasi- Id gapless MSL can be discussed as in the spinless case 
(see also Appendix [All. The space-time gauge field fluctua- 
tions are suppressed by the destructive interference arising 
from the incommensurate momenta carried by the fermion 
fields. Thus, the system retains three gapless modes, but the 
local energy observable obtains new oscillating contributions. 

We can also consider directly allowed perturbations going 
beyond the exactly solvable model. For example, rf terms 
on the vertical links {ij) can be expressed as a product of 
three c-fermion strings, one for each flavor, and will oscillate 
at wavevector 3kp with power law X^^^^ in the free parton 
case. This is consistent with the schematic analysis in Ap- 
pendix |A] extended to multiple parton fields, where a vison 
can be seen as introducing a half-vortex for each flavor. The 
described low-energy theory is hence stable to generic per- 
turbations in the sense of retaining the gapless fields, while 
the local energy observable that is symmetric under the leg- 
interchange obtains additional contributions oscillating at 3k p 
(which in turn induces new contributions to other observables 
as discussed in AppendixlBl). 



IV. DISCUSSION 

Motivated by recent proposal of SU(2)-invariant Majorana 
Spin Liquids by Biswas et ali^ and the realization of the 
SU(2) MSL in an exactly solvable modelr^i^^^ we studied 
the MOL and SU(2) MSL on the 2-leg ladder Perturbing 
away from the exactly solvable points, in the MOL, there is 
only a strictly marginal four-fermion interaction and hence it 
is stable to residual interactions. In the SU(2) MSL, there 
are several allowed four-fermion terms, but it is stable against 
these in a large parameter regime. Furthermore, we also show 
that such gapless Majorana liquids persist against Z2 gauge 
field fluctuations. Some time ago, Shastry and Ser^ stud- 
ied an SU(2) MSL for a Id Heisenberg chain at mean field 
level. Our description of the microscopically realized quasi- 
Id SU(2) MSL can be viewed as providing a theory beyond 
mean field for more general such states and distinguishes them 
from the Bethe phase of the Id Heisenberg chain. The stable 
MOL and SU(2) MSL phases that we find are new quasi-ld 
phases, and we suggest numerical studies such as Density Ma- 
trix Renormalization Group (DMRG)^ to test our theoretical 
ideas of their stability. The DMRG studies can also determine 
the Luttinger parameters of the fixed-point MOL and SU(2) 
MSL theories. 

The presence of gapless matter fields is the key against con- 
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fining effects of Z2 gauge field fluctuations in (1+1)D, see 
Appendix lAl Without such gapless matter, the gapped phases 
realized in Kitaev-type models on 2-leg ladders in our model 
are likely unstable to general generic perturbations, and this 
prediction can be checked by DMRG studies. This is remi- 
niscent of a picture where gapless matter fields can suppress 
monopoles in a (2+l)D compact electrodynamics and thus 
make gapless U(l) spin liquids with sufficiently many Dirac 
points or with Fermi surfaces stable j'-^-'° while gapped U(l) 
spin liquids would be unstable to confinement in (2+l)D. An 
interesting finding is that allowing Z-i gauge fluctuations in 
our quasi-ld Majorana liquids leads to new contributions to 
various observables, with different characteristic wavevectors 
and potentially slower power laws compared to the mean field, 
cf. Appendix IB] 

Let us discuss possible extensions of this work. Through- 
out, we focused on the MSL phase in which all couplings 
of the residual interactions, Eq. ( |37] |. converge to finite fixed 
point values in RG thinking. In principle, one can analyze 
situations where some of the residual interactions are rele- 
vant and explore possible nearby phases and characterize their 
properties using the observables listed in Appendix |B] Such 
theoretical analysis combined with DMRG studies^ can give 
a complete phase diagram. 

As discussed in Biswas et alM- and in our earlier work^^ 
the effects of Zeeman field on the SU(2) MSL are interest- 
ing. The Zeeman magnetic field only couples to and 
fermions, and we can define /| = (/^t ± ip^)/^ which 
carry 5^ = ±1, while carries S*^ = and remains un- 
altered. In the presence of the Zeeman field, the spin SU(2) 
rotation symmetry is broken and only 5^ is conserved. In 
Appendix ICl we write down general four-fermion interactions 
based on symmetry arguments and perform weak coupling 
RG analysis. Our RG equations ( IC3b - (IC7b interestingly show 
that instabilities only occur in the channel but not in the 
channel. Hence, the partons are always gapless no matter 
how large the field is and can give metal-like contribution to 
specific heat and thermal conductivity, which is qualitatively 
similar to what we found previously in our 2d MSL modeli^ 

Last but not least, it is intriguing to understand how the 
ladder descendants of the MOL and SU(2) MSL relate to the 
mother 2d phases. A systematic way to access these could 
be via increasing the number of legs. It seems difficult to in- 
crease the number of legs in our toy 2-leg square ladder model 
while maintaining the spin SU(2) symmetry of the MSL, but 
actually it can be achieved if we consider decorated square 
ladder,--'-"' One more interesting direction is to consider new 
types of SU(2)-invariant spin liquid wave functions motivated 
by the Kitaev-like SU(2) MSL writing of the spin operators 
and search for more reaUstic models in Id and 2d that may 
harbor such states. 



Appendix A: Stability of gapless U(l) matter against Z2 gauge 
field fluctuations in (1+1)D 

We need to address the issue whether the gapless parton 
field picture is stable against allowing Z2 gauge field fluctu- 



ations. It is well-known that the simplest so-called even Z2 
gauge theory is confining in (1h-1)D; this persists also in the 
presence of gapped matter fields, and quasi-ld Kitaev-type 
models with gapped partons would suffer from this instabil- 
ity. We will argue, however, that gapless parton fields can 
eliminate this instability, particularly when they carry incom- 
mensurate momenta. 

We first give a heuristic argument. Let us consider the sim- 
plest model of a Z2 gauge field coupled to a U( 1 ) matter field, 
with (1h-1)D action 

5 = (7jk cos - (/)fc) - cri2Cr230-340'41-(Al) 

For K 00, we choose the gauge ajk = 1 and obtain an 
XY model in the (p variables. There is a Kosterlitz-Thouless 
transition at some critical /3c and gapless phase for (3 > /3c. 
Now, let us consider large K and large /3 limit. Starting 
with no Z2 fluxes and no vortices, since both cr and <f> are al- 
most fixed, the insertion of a Z2 flux ("vison") can be treated 
as creating a 7r-vortex in the (f>. Explicitly, we can rewrite 
(Tjk COS - (?!)fe) = cos - (?!)fe - 7r(l - crjfc)/2]. The vi- 
son insertion can be carried out by changing ajk from 1 to 
— 1 on a cut from infinity to the vison location. This is a tt- 
phase cut for the cj) variables and can be best accomodated by 
a gradual winding by tt as we go around the vison from one 
side of the cut to the other; hence, we get a half-vortex in the 
<j>. We expect that for sufficiently large /?, the half-vortex in- 
sertions are irrelevant because of their high energy cost, which 
means we have a phase without proliferation of half-vortices, 
and then we do not need to worry about the dynamics of the 
Z2 gauge field which could potentially produce confinement. 

Thus, it is possible to avoid confinement of (1h-1)D Z2 
gauge fields if we have gapless matter field. For several gap- 
less matter fields, there is a proportional increase in the en- 
ergy cost of the vison insertion and hence its irrelevance. The 
above argument is valid for matter fields at integer filling. It is 
well-known that vortices in (1h-1)D U(1) systems can be fur- 
ther suppressed if the matter field is at non-integer filling due 
to Berry phase effects, and such a suppression is complete for 
incommensurate matter density. Heuristically, we expect the 
vison insertions to obtain similar Berry phases as half-vortices 
and hence to also experience complete suppression at incom- 
mensurate density. We present a more formal derivationii tai- 
lored to our needs below. 

We consider a general Z2 gauge theory plus C/(l) matter 
field (represented by quantum rotors) on a d-dimensional cu- 
bic lattice with a Hamiltonian" 

H = —t ^ a^j., cos {(f)r — 0r') + ^ 5^^"'' ~ "■'^ 
{rr') r 

-K J2 - r ^ a^„, , (A2) 

□ {rr') 

where fir is the number operator conjugate to the phase 0^ at 
site r and n is the average density. The Hilbert space con- 
straint is 

e""- \[ a^r' = 1 • (A3) 
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We proceed to treat the system using standard Euclidean path After standard development of the path integral for the Ising 
integral formalism in the cr^-(/) basis. We implement the con- gauge fields, we can write the partition function as 
straint at each site r and temporal coordinate t by using the 
identity 



1 



E 



A(r,T)=±l 



J 



J2 r-DMr) E 

_N _M "'0 r__ 



{Sf ,(T);A(r,r)} 



xe 2 



(A4) 



Here we used S^.^., to denote eigenvalues of a^^., on the spatial over all spatial and temporal plackets, Kp = {Xspat, Kt}, 
links and elevated the auxiliary fields A(7', t) to become Ising with Kspat = K5t and tanhX^ = e~'^^^'^ . 
gauge fields on the temporal links, 5"^^^ ^y^^ t+St) = '^(^' ^) Now we can use a variant of XY duality transformation^— 

(we use either field notation where more convenient); J^p is 8° *^ ^"'^ " variables to real-valued "currents" 

jspat = {jr,r+eiJr,r+e2,---,.ir,r+ea} (where efc=i...d rep- 
resent unit lattice vectors) and jr appearing as follows: 



^tirS^ ,(r)cos[0,(r)-0^,(T)] 



+25 -i|^k,(r)-0,-(r)+,r 

Prr' (t) — — OO 

+00 ^+oo 



E 

p^^/ (r) — — oo 



djrr'{T)e 



^fe^ + y..' (^) 0,.' (r)-0,(T)+7r 



l-S; , (t 



-27rp,,,(T) 



^ F[nr[T)]= dMr,r) ^ e-^^--('-'-)-2-P^('^'^)F[jV(r, r)] 

/ \ — OO / \ 



(A5) 
(A6) 



'n.r{T) — — 0O 



Pt {r,r) — — cxd 



r 



In the first Une, we approximated the left hand side by a stan- where p[^^ is one representative of a class; the results do not 

dard Villain form; we also dropped constant numerical fac- depend on the specific choices of but only on the "vortic- 

tors throughout. For short-hand, we write space-time points jjies" ^ y^p^ _ y^p^ characterizing the classes. Us- 

asj = (r, r) andjlefine space-time vector = jjjg ^jje Ni variables, we can extend the 0j integrations to 

{Pspat, Pt}i with Pspat = {Pr,r+en Pr,r+e2 t ■ ■ ■ t Pr.r+Cd} ■ ( — OO, +Oo) and obtain 

Then we can divide configurations {pi^i} into classes Cp 
equivalent under integer- valued gauge transformations pi^ — > 
Pifi + '^fj.Ni and perform the configuration summation as 



+ 00 



{Pin} = -c 



/OO 
'D]spa.t'Djr (5(V • jspat + Vrjr = 0) X e 
_ , . . , -OO 



Ep -^P*5'i2'5'23'^34'^41 



{S; ,(r);A(r,r)} Cp 



'27rp<°',(T) 



xe 

Xe-^x,,-^[j4'-.r)-n]"+»E.,,.jW^.r)[7r ^-^f-"' -27rpt")(r,r)] 
I 
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(A7) 



The above result holds in general (dH-l)Dri and from now on we specialize to (1h-1)D system. We solve the current con- 
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servation condition by writing = n^^^ = ^^^^i^, with 
^(x, r) = -Kfix, X being the spatial coordinate on the dual lat- 
tice, and = = The dual field e encodes 
coarse-grained fluctuations in the particle number. 



We have only temporal plackets, on which we define "vor- 
ticity" q = W X p ^ V xPt — ^ tPx and "vison number" 
j^vison ^ V X (1 - 5^)/2 mod 2 = or 1 corresponding to 
— 1 or —1. We can absorb any modulo 2 shifts 
from 71^'*^°" by redefining q and write the partition function as 



This is the main result, which we can now analyze in a num- 
ber of standard ways. We can integrate out the field 6 and ob- 
tain a Coulomb gas representation. In the absence of the Z-2 
gauge field (e.g., X — > oo and 71^'^°" = 0), we get famil- 
iar integer- valued charges q representing vortices of the U{\) 
matter system. On the other hand, for any finite K we get 
effectively half-integer charges m ~ q — ^n^^^°''^ G i x Z 
with only short-scale energetics difference between integer 
and half-integer charges. We also see Berry phases e'^Trna- 
for a vortex insertion in the presence of non-zero background 
density and halving of the Berry phase for a vison insertion. 
Alternatively, we can consider postulating some local energet- 
ics penalty for large values of m and perfom the summation 
over to obtain terms like 



Ai/2 cos(6' + nnx) + Ai cos(26' + 2'Knx) + 



(A9) 



where we ommitted possible phase shifts in the cosines for 
brevity. The Ai term is the familiar term in the dual sine- 
Gordon theory for a Luttinger liquid of bosons that represents 
allowing vortices, while the A1/2 term can be now interpreted 
as effectively allowing half-vortices if the matter is coupled 
to Z2 gauge fields. Crucially, both vortices and visons expe- 
rience destructive interference effects for incommensurate n. 
On the other hand, for commensurate n the vison insertions 
can still be rendered irrelevant by going deep enough into the 
Luttinger phase or increasing the number of gapless fields as 
discussed below. 

We can generalize the above result to the case with sev- 
eral matter fields coupled to the same Z2 gauge field 
by replacing the Berry phase 2{9 + T:nx) x {q — iri^"'°'^) 
wifli J2a + T^fiax) X (g„ - iri'''"°'^). Here the sum- 
mation over vison numbers leads effectively to terms like 
Ai/2 cos(X]q + Tr^^nax). We can see that for three 
identical flavors with incommensurate n as happens in the 
SU(2)-invariant MSL, the destructive interference effects will 
wash out any vison insertions (including any combinations 
with non- vison terms). 

Looking back at the one-component case, we could ratio- 
nalize the above structure more quickly by thinking about the 
theory Eq. (IA2b as coming from a formal splitting of some 
physical boson field e*'^p''*'=' into two halves:" schematically, 
gi^phys _ gi20 xjjen the described gapless phase can be 
thought of as a (1h-1)D analogue of the "Higgs phase" that is 
expectedii to reproduce the conventional "superfluid" (here. 



quasi-long-range ordered) phase of the physical bosons. In- 
deed, in the derived harmonic liquid description in terms of 
the dual field 6, we can change to new variable (?phys — 0/2 
canonically dual to 0phys and note that the identified vison in- 
sertion operator e'^ = e'^^p''''^ is the same as the conventional 
vortex insertion in (/)phys ■ We still like to show the above more 
formal derivation as it is not tied to the specific origin of the 
parton field cj). For example, in Sec. Ill AI the conjugate pair 
{0, 6} arose from bosonizing the long-wavelength fermionic 
parton Hamiltonian, and we can continue using these fields in 
calculations but remember to include the Z2 gauge fluctuation 
effects by allowing local energy terms like A1/2 cos(6' + 7rnx). 
The same formal treatment also holds transparently for the 
multi-flavor generalization where the parton fields provide 
a very convenient description of the unconventional gapless 
phase, which has the same number of gapless modes as in the 
parton mean-field, but with the identified new contributions to 
the local energy once we go beyond the mean field and include 
Z2 gauge field fluctuations. 



Appendix B: Fixed-point theory and observables in the SU(2) 
Majorana spin liquid 

We use Bosonization to re-express the low energy fermion 
operators. 



.1 P — 'lad , 

with canonical conjugate boson fields: 

[^^{x),^p{x')] = [e^{x),epix')] = o, 

[(Pa{x),6j3{x')] ^ inSapQix - x'), 



(Bl) 

(B2) 
(B3) 



where Q{x) is the Heaviside step function and we have intro- 
duced Klein factors, the Majorana fermions with {ija, ?/^} = 
2(5q/3, which assure that the fermion fields with different fla- 
vors anti-commute with one another 

According to the RG analysis in Sec. [nil at the fixed point 
of the stable SU(2) MSL phase, only the coupling Up is strictly 
marginal and will renormalize the Luttinger parameter g in the 
"charge" sector. The effective bosonized Lagrangian is 



5(7(2) 
MSL 



1 



Vn 



Vp{d, 



E 

p,= l,2 



1 

2^ 



-{dre^f + v{dxep 



(B4) 
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where we defined 



symmetric (a) under the leg interchange: 



1 

7!' 
1 

7!' 
1 

7!' 



-Oy), 

+ 6y- 26,) , 



(B5) 
(B6) 
(B7) 



and similarly for the ip-s, which preserves the commutation 
relations, Eqs. ( IB2b - (IB3b . Stability against the W4 term in 
Eq. (|37] | requires g < 1- 

For the observables characterizing the SU(2) MSL phase, 
as discussed in Ref. IIH we can use spin operators. 



S - ^ 



bond energy operators. 



and spin-nematic operators 



(B8) 



(B9) 



(BIO) 



P=R/L 

-*E^"'V^/2, 

~* E ^R^L^ 
- ~* E ^R^L^ 



kFR+kpL 



na 

kpR+kpL 



kp 

QCt.S 



FR+kpL 



kpR—^FL 



/3,7 



■7 
R^ 



Q 



FR — i^FL 
''kFR — kFL 



E f^L f \ 







E ■I'' 

Jl Jr. "I" Jl JRy 



■'2kf 
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(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 

(B18) 

(B19) 
(B20) 

(B21) 

(B22) 

(B23) 
(B24) 



The latter can be related to the usual traceless rank two 
quadrupolar tensor defined as 



■]k 



2 \'^J ^k 



Qp oa 
'^k 



-5''f'{S,-Sk), (BID 



through P+ = Q^j^^Q 



Jk 



^^Q?k- 



We expand the observables in terms of the continuum com- 
plex fermion fields and organize according to the momentum 
and the leg interchange symmetry, i.e. symmetric (s) or anti- 



with 5"q = S'q , etc., and C*/" observables mean sym- 
metric or anti-symmetric under the leg interchange. If the 
TRS is broken explicitly as in Appendix |D] all the above 
momenta are distinct. With TRS, kpL = —kpp, we have 
coincident momenta kpR + kpL = and kpp — kpp = 
^kpp = —2kpL. Strictly speaking, with TRS, we should 
define 0^^„ = '^fcj.H+fcFL + H. c, ins tead of Eqs. ( IBT6] i- 
( |B19t ; similarly, instead of Eq. ( IB24b . we should define 

'^2fe" ~ '^2kFR ~^ ^-tkFL ■ '■^^ present case, the listed terms 
with such equal momenta transform differently under leg in- 
terchange, which is encoded in the above definitions. 
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The bosonized forms at Q = are: 



Sq=o = cos ^= j cos 

- 4^n n cos ( + cos (^1^1^ 

'^Q=0 — ^^VxVz COS ^ ^ J COS ^ ^ 



5^ 



4M/j,r/x cos (■\/2(pi) cos 
V3 



f V3^2 - "yfl ^ ( V^02 - 6»i , 

,(B26) 
(B27) 

(B28) 
(B29) 
(B30) 

(B31) 
(B32) 

,(B33) 

The corresponding scaUng dimension in the fixed-point theory 
Eq. dm is 

A[S^Q^,] = A[6^^o] = mi^a = 1, (B35) 

which is not modified by the strictly marginal interactions. 
The bosonized forms at = kpn + kpL are; 



1 [2, 



Qxx,s 
Q=0 



Sq=o = ^'i-VyVx COS (^V20i^ sin (^%/2v3i 
S^to = 4i77^% cos ( ^^^ sin ( ^'^^ '^^ 
^0=0 = '*'^VzVx cos I I sm 



2^VzVy<^^ ^ v^'cOS — 

V v2 
fV302 + Oi 



S^n" = 2i77,r,,e^(7f'^''-^+^) cos 



\ 72 



' cos(%/26'i), 



2e»\/f¥'2 cos(y2^i) + e"*^^'^^ 



s: 



s: 



Ql^^^- = 2rjxVze^^^'''-^-'^^' sin 
The corresponding scaling dimension is 



V V2 
V2 



2 1 

3 ^ 



(B36) 

(B37) 

(B38) 
(B39) 

(B40) 

(B41) 

(B42) 
(B43) 

(B44) 
(B45) 

(B46) 



The bosonized forms at Q- = kpR — kpL are: 

S^l = 2^,.,,e'(^^-^-^) cos ( ^^1^ 



V2 



V V2 



Sq_ = 2irjy'qxe'"- ^ 



cos(^/2v3l), 



= le ^/3 



2eV3S2 cos(y26li) + e 
gi(y|e2-y2ei) _|_g-*2y|( 



42^462 



2ze'^7f '^p-v 3-2; cos(V26li), 



qzz.s _ 
■\xy.s 



Qyz.s o- 
^_ = 2iri^r]ye 



*'■^/3*'P ^/6 v^-* gin ' ^ ^-i- 



V x/2 



i(2g 83^_^lL) , /\/3(y52 + Vl 



The corresponding scaling dimension is 



(B47) 

(B48) 

(B49) 
(B50) 

(B51) 

(B52) 

(B53) 
(B54) 

,(B55) 
.(B56) 

(B57) 



The bosonized forms at the 2k pp are: 

S'2fc" = 2i?7^?/ye'^^'^''"75-^)e'-^(^^p-;;%-^\ (B58) 



S2kpp = 2i?7y?7^e'^7S 
where P = R/L = ±. 



« ( TT? p + 2 ) 4 P ( -f, e p + e 2 ) 



1 

\ (B60) 



(B61) 



1 5 J_ 
3^3^ 3.g' 

We can see that when g ~ 1, each scaling dimension is 1, 
the value in the exactly solvable models with non-interacting 
partons. In the stable SU(2) MSL, we require g <1 and hence 

A[OqJ < A[Oq^o] < A[02fc,,] < A[OqJ (B62) 

Besides the observables constructed out of local fermion fields 
discussed above, there are local physical observables that re- 
quire non-local expressions in terms of fermion fields similar 
to the string operator defined in Eq. ( |28] l. In this SU(2) case, 
we can consider the "rung energy" operator which is symmet- 
ric under leg interchange, 

e{X) = T%X,l)T%XA). (B63) 

Considering correlation function of such an operator similar 
to Eq. ( |26] | in the spinless case, we can write schematically in 
our gauge 

r^(X,l)r^(X,4) - 

- n n("l)^"(^''l)^"(^''4)"="(^''2)c"(X',3).(B64) 

X'<X a 
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Such non-local operator in fermionic language seems very in- 
tractable but the expression can be greatly simplified under 
Bosonization, 



,umt cell 



n n ^"(^'' 4)c"(X', 2)c"(X', 3) 

X'<X a 

^ p±iT.JS{X)+7rn^X] ^ ±^[^/30p+3fcFX] 



(B65) 



where we used the definition of 9p in Eq. iB5i . kpR = kp 
and fia = kp/TT is the average density of a-species fermion. 
Thus, we can write a contribution to the leg-symmetric energy 
observable as 



<^3kf 



(B66) 



with scaHng dimension A[e3fcp„] = ^ and 63^^^ = 
Egj, . We can also consider other rung energy operator 
such as r*'(X, l)Ty(X, 4), but the long-wavelength descrip- 
tion of such an operator is qualitatively the same as the 
above r^(X, l)r^(X, 4). Finally, these local energy observ- 
ables can be combined with any observables listed earlier to 
produce further critical operators with potentially enhanced 
scaUng dimension, e.g. Ofc^^+afc^z, ^■ikpi.OQ_ with 

^\0kFR^2kF,\ = 3 



- I + ^ and S^f.^^^^f.^^ 



<^3kpLS2ki. 



_a_ 

12 



j_ 
3g- 



Appendix C: Zeeman magnetic field effects on tlie SU(2) 
Majorana spin liquid 

In the SU(2) MSL phase, Zeeman magnetic field only 
couples to and fermions, and we can define f]. = 
(/^■f ± i/^^)/\/2 which carry ~ ±1 and get Zeeman- 
shifted, while f^^ carries 5*^ = that remains unaltered. 
The spin SU(2) rotation symmetry is broken and only is 
conserved. Using symmetry arguments, we can write general 
four-fermion perturbations in terms of long-wavelength right- 
moving and left-moving complex fermions as 



-W^^{.f+,Rf+,Lf- 



differential RG equations are 



with pf,,p = f^^f^, p. 



A^ 



w 



+- 



A+ 













{w+-f 






? 




{w+ 


-? 




tt{v+ - 


\-V- 






■A+ 


+ 


2n 


V- 




' = X~ 


z _ 


= 0. 



(CD 

,,j/„,L + H.c.) , (C2) 
z, and P = R/L. The 

(C3) 
(C4) 
(C5) 



+ 



4A+- 



(C6) 
(C7) 



h:: 



hz 



X 3 ^ 



h::i 



h z' 



1 



(a) 
6* 



1 



X 2 y 







2 


Band-1\ 






1 


Band-2 y 


— K 






1 ^/ 



(b) 



FIG. 6. (a) Graphical representation of the exactly solvable Kitaev- 
type model with time reversal breaking (TRB) introduced by hand 
and its solution in the zero flux sector, (b) Complex fermion 
spectrum, Eq. for the Majorana spin liquid with TRB with 
{J., Jy, J., J^, h} = {1.2, 0.8, 1.0, 1.1, 0.5}. 



Here O = dO / d£, where £ is logarithm of the length scale 
and v± represent Fermi velocities of the bands. We see 
that the MSL is stable if 



A^ 



A- 



4A^ 



> 0. 



(C8) 



Comparing the RG equations (IC3b -( IC7b in the presence of the 
Zeeman magnetic field with those Eqs. (I4ni-(l44li without the 
Zeeman field, we see that the instabilities in the "spin" sec- 
tor, Ua-i and Ua-2, are removed by the magnetic field, and the 
couplings that contain both and do not flow (the reason 
is that interactions that could cause these to flow do not con- 
serve S*^ and thus are not allowed). An interesting fact about 
these RG equations is that the instabilities only occur in the 

fermion but not in the channel. Hence, the gapless 
partons are always gapless no matter how large the Zeeman 
magnetic field is and always give metal-hke contribution to 
specific heat and thermal conductivity. 



Appendix D: SU(2) Majorana spin liquid with Time Reversal 
Breaking (TRB) 

In this Appendix, we will break the time reversal symmetry 
explicitly by including a term. 



T~iTRB 



h 



E 



(TMr|-r^r|rr)(a3-ai) (Dl) 



+ (^1^1 ^4 - rlTlT^) ((74 • a2) .(D2) 
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Later we will see that such terms reduce the number of four- 
fermion interactions due to momentum conservation. Using 
the Majorana representation, this term can be rephrased as 



(U34U4I 



ui2U2:i) ^3C? (D3) 



(U41U12 + U23"34) ^ C4 C2 



(D4) 



The graphical representation is shown in Fig. |6(a)| Before 
we proceed, we remark that in this case with TRB, we do 
not need any symmetry to protect the gaplessness, unlike the 
time reversal invariant case. The bilinear term IftL that could 
open a gap is not allowed in the Hamiltonian due to momen- 
tum conservation, see below. For illustration and simplicity, 
we proceed to take the same parameters as in Sec.|III]and in- 
clude /i, {J^, Jy, J,, j;, h] = {1.2, 0.8, 1.0, 1.1, 0.5}. 
The complex fermion spectrum is shown in Fig. |6(b)[ and we 
can clearly see that due to the presence of the Time-Reversal 
Breaking term, there is no Right-Left symmetry anymore (i.e. 
kpL 7^ —kpji). In the weak-coupling regime, the general 
four-fermion interactions can be written as 



K 



TRB 



UhJrJl — UctiJr ■ Jl 



Ua2l\i^lR.L, (D5) 



where jTp, jTp, and Ij^l are defined in Eqs. (l38Tl-(l40li. We 
can see that the number of allowed interactions is reduced be- 
cause there is no special relation between kpR and fc^L and 
additional terms are forbidden by momentum conservation. 

The weak-coupling differential RG equations in this case 



, _ 1 
, 1 



[?i^2 + 2iio-lWa-2] , 



(D6) 
(D7) 
(D8) 



We can give a qualitative description of the stable flows3 
If Ua\ > and Ucri + > 0, the coupHngs Uai,2 are 
marginally irrelevant and flow to zero, it*]^ = u*2 = 0. 
The coupling tip approaches a fixed value, u*, and is strictly 
marginal; unlike the time reversal symmetric case in Sec.HIIl 
there is no condition on the sign of u*. We conclude that 
the SU(2) MSL with explicit time reversal breaking is sta- 
ble in a wide regime of parameters. We also note that, even 
though initially there is no conservation of the /-fermions in 
this model, breaking TRS leads iokpL 7^ —kpR and prohibits 
four-fermion interactions such as f"f^f^f^ and f^^f^f^f^, 
so the fermion conservation emerges at low energy. We note 
that if we rewrite the couplings as 



Ua2 



7r(uH + Vl) 



2V2 



-91, 



2V2 

the RG equations can be rephrased as 



(51+32), 



9p 
91 

92 



3Up + Ua2 = 0, 



1 



(3.9? + 2.91,92) , 
2.91.92) 



2V2 



(D9) 
(DIG) 

(DID 
(D12) 

(D13) 



The last two equation are exactly the same as one-loop RG 
equations in an SU(3) WZW model in Ref. |M Note that in 
the SU(2) MSL the "charge" (p) sector also remains gapless, 
cf. Appendix IB] 
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